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Abstract 

We analyse the chaotic motion and its shape dependence in a piecewise linear map 
using Fujisaka's characteristic function method. The map is a generalization of the one 
introduced by R. Artuso. Exact expressions for diffusion coefficient are obtained giving 
previously obtained results as special cases. Fluctuation spectrum relating to probability 
density function is obtained in a parametric form. We also give limiting forms of the 
above quantities. Dependence of diffusion coefficient and probability density function on 
the shape of the map is examined. 
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Deterministic diffusion is a well known phenomenon in spatially extended 1-D maps 
[1-8]. It has been proposed as a possible mechanism to account for the behaviour of 
Josephson junctions in the presence of microwave radiation [9] and of parametrically 
driven oscillators [10]. In Hamiltonian dynamics [11] also, transport due to chaos is sig- 
nificant because of its applications in celestial mechanics, confinement problems and so 
on. Recently, some exactly solvable models have been analysed [6-8]. The only aim in 
these studies is the evaluation of the exact diffusion coefficient using a cycle expansion 
technique [12]. It is a well known fact that the chaotic dynamics in spatially extended 
maps has two complementary aspects - diffusion and intermittency. These are related to 
the probability distribution which is approximately Gaussian by central limit theorem. 
Fujisaka's characteristic function method is a useful tool for analyzing both these aspects 
of stochasticity in such maps. In this brief report, we apply the characteristic function 
formalism [13-14] to analyse the chaotic motion in a generalized piecewise linear (GPL) 
map with a variable shape. It is a generalization of the exactly solvable model in ref. [6] 
allowing analytical study. Exact expression for diffusion coefficient and a parametric 
representation for the fluctuation spectrum relating to the probability density function 
(PDF) are obtained. Generalization permits the study of the dependence of these quanti- 
ties on the shape of the map. We notice that GPL map with flat peaks is more suited to 
describe systems exhibiting intermittency in time. The generalization brings the map in 
ref. [6] nearer to sinusoidal maps studied numerically in ref. [1]. A similar shape dependent 
piecewise linear model has been examined in ref. [2] from the point of view of correlation 
times. 

Chaos-induced diffusion systems have a general form [4,13] 

X t+1 = X t + P r (X t ) = Y r (X t ), P r (X + 1) = P r (X) (1) 

where r is a control parameter. The sinusoidal map P r (X) = rsin(2HX) is an example 
[1]. After the decomposition X t = N t + x t where N t is the cell number measured in which 
X t is located and x t (0 < x t < 1) the distance measured from the relative origin X = N t , 
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eq.(l) can be uniquely rewritten as two dynamical laws: 



t+i — 



N t + A(x t ), x t+1 = f(x t ) 



(2) 



Here A(x) is the jumping number defined as the largest integer smaller than x + P r (x) 
and is free from N t and f(x) = x + P r (x) — A(x), satisfying < f(x) < 1. f(x) is the 
reduced map of the extended map (1). 

We analyse a piecewise linear map with variable shape of the type in fig.(l). In the gen- 
eral case, the map consists of linear segments with slopes ±TOj, i — 0, 1, • • • , h, rrii < m^. 
For the cells on the bisector, the slope magnitude is m . For the i th cell above and below 
this cell on bisector, the slope magnitude changes to m 8 . The reduced map consists of k — 
Ah + 3 linear segments. For k increasing from 1 to Ah + 3, these line segments have slopes 

mo, mi, m 2 , ■ • • m h , -m h , -m h - U • • • 0, -mi, -ra 2 , m h , m h , m h _i, • • •m 2 ,m 1 ,m , and 

mjS satisfy the relation 



The extended map can be generated from the reduced map by giving suitable jump 
numbers A(x) (constant for a line segment). These (from left) are 0, 1, 2, • • • h, h, h — 

l,---2,l,0,-l,-2,---,-/i,-/i,-(/i- 1), 2,-1,0. Figs.(l) and (2) show the map 

and the reduced map for h — 1. 

The map (1) can be studied using the characteristic function formalism [13-14]. In 
this, the dynamics of A t governed by A t+ \ = B(x t )A t (t = 0, 1,2, •••) with A = 1 is 
studied. B(x t ) is a steady function of x t which evolves according to the chaotic map 
%t+i = f(x t )(0 < x t < 1) [13]. Equivalently, one can consider the dynamics of the local 
time average of a time series a t = \ Y?j=i l n B{ x j) [14]. Map (1) can be treated by putting 
A t = exp{N t - N );B{x) = exp(A(x)). We put N = 0. Then a t = The long time 
dynamics of N t can be studied using Fujisaka's characteristic function 
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(3) 



A, = - lim -ln[{ex V (qN t ))} 

q t— »oo i 



(4) 
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(exp(qN t )) is the average over a steady ensemble and is the q-order moment of exp(iV t ). 
One can expand \ q in the series of cumulants. The expansion converges for \q\ < c, c 
being the convergence radius. In this case, \ q can be approximated as 

\ q = A + D q (5) 
Ao = «oo = hm — 

t^oo t 

where Ao is the drift velocity. D is the diffusion coefficient given by 

a t =((N t -X t) 2 )^2Dt (6) 

for large values of t. a t is the variance of N t . Asymptotic PDF of a t has a Gaussian 
component (central limit theorem) and a non-Gaussian component. For \q\ << c, the 
moment < exp(g N t ) > is determined by the Gaussian component (diffusion branch of q) 
and for \q\ » c, it is determined by the non-Gaussian component (intermittency branch 
of q). The PDF pt(a) that a t takes values between a and a + da can be obtained as 

p t (a) ~ exp[-a(a)t] (7) 

a(a) being the fluctuation spectrum, pt{a) — > 5(a — a^) as t — > oo. pt(«) can be obtained 
from Ag in parametric form using the Legendre transform 

^(«) = ? 2 ^^ (8) 

We first consider the case with /i = 1. The reduced map consists of 7 line segments 
with slopes (from left) m , mi, —mi, —m , —mi,mi, m . These satisfy (3). The Frobenius- 
Perron operator H is defined by 

where yt is the k th solution of f(yk) = x and f'(x) = ^f{x). m k is the slope of the k th 
line segment of the reduced map. From eq.(9) we note that the invariant density p*(x) is 



uniform (p*(x) = 1) in the interval < x < 1 [Hp*(x) = p*(x)]. The Lyapunov exponent 
A can be obtained as 

3 4 

A = — ln(m ) H ln(rai) (10) 

m mi 

Since m 8 > 1, we note that A > and therefore the reduced map is always chaotic. 
Characteristic function \ q can be evaluated [13] using the linear operator defined by Mori 
et al [15] 

UF(x) = -±^H\p*(x)F(x)] (11) 

A = - lim - ln( e qAx H e qAx H e qAx H e qAx ) (12) 
q q t-oo t x ' . ' 1 v 1 

For our model H = H. A(x) which are constant over a line segment are (from left) 
0, +1, +1, 0, -1, -1, 0. Hence we get from Eqs.(9), (11) and (12) 

H e qAx = — + — e q + —e- q 
m mi mi 

= — + — cosh(g) (13) 
m mi 

\ q = -ln[— + — cosh(g)] (14) 
q m mi 

The result can be generalized for integer values of h. Again, slopes satisfy relation (3). 
Frobenius- Perron operator again leads to uniform invariant density p*(x) = 1. Lyapunov 
exponent A and characteristic function \ q are given by 

3 h 4 

A = — ln(m ) + V — ]n(mi) (15) 
mo ~i m i 

A = - ln[— + Y — cosh(i g)l (16) 
g m ^ mi 

Map is fully chaotic since A > 0. The drift velocity A = 0, always. The diffusion 
coefficient D is 



D = lim Aa, = V — (17) 

The fluctuation spectrum a(a) can be got in the parametric form using Eq.(8) 

Etx^sinh^g) 



a 



— + £?-i — cosh(ig) 

ran 1 m, V M/ 



mo i mi 

Eti^anMi?) . , 3 » 4 



^) = «X^X^ -M- + g-cosh( i9) ] (18) 



mo 



i=l 



q = gives a = 0; <r(a) = 0. If +g gives +a, —q will give —a without changing a (a). It 
can also be noted that maximum value of a is obtained by putting q — > oo. We have 

aw = h a(a max ) = ln(^y-) (19) 

In the special case when all m^s are equal (=m ), eq.(17) can be summed to obtain a 
closed form expression for D. In this case, eq.(3) gives m = 3 + Ah. 

_ fc(fe + l)(2fr+l) 
^" 3(4/^ + 3) 1 Uj 

With /i + 1 = (3, 

((3 - 1)/3(2P - 1) 

3(4/5-1) 1 j 

which reduces to -D = 2/7 for /3 = 2 giving results obtained previously in ref.(6). 

Closed form expression can be got also for the special case mj/mj-i = r = a constant, 
< r < 1. Then rrii = m^r 1 . From eq.(3) 

4(1 ~r h ) , . 

r ft (l — r) 

For every r between and 1, the model becomes an exactly solvable case. From eq.(17) 

n _2[/i 2 + (1 - 2h - 2h 2 )r + (h + 1) V - r h+1 - r h+2 ] 

~ [3r ft (l-r)+4(l-r ft )](l-r) 2 (23) 

Limiting forms of the above quantities can be obtained for a constant h as the peak 
shape becomes maximum flat. These can be arrived at by taking limit r — > 0. From 
eq.(23), D behaves like 



h 2 

D = \ (24) 

The above limit can also be obtained by putting — > oo(? = 0, • • • h — 1) and m ft — > 4. 
Applying this we get the following limits 

lim \ q = - ln[cosh(g/i)] (25) 
lima = /itanh(g/i) (26) 

lima(a) = ln[^-^]V 2[ ^]^ (27) 

In fig. (3) we plot diffusion coefficient versus r for /i = 2. It can be observed that D 
increases with increasing flatness of the peak shape. D varies from 0.9 to 2 when r is 
varied from 1 to 0. Increasing h, keeping rrii = m (? = 1, 2 • • • h) appears to have more 
influence on increasing D. This is because _D = 2/7 = 0.29for/i = l, whereas it goes to 
0.909 for h = 3 (eq.(20). 

The probability distribution function for N t , the distance from the origin, can be 
obtained using the fluctuation spectrum a (a). From eq.(7), we have 

1 N 

^(jV)~-exp[-<r(-)f] (28) 

p t (N) being the PDF that N t takes values between iV and iV + dN. This PDF is ap- 
proximately Gaussian by central limit theorem. In the exactly normal case \ q is given by 
eq.(5) and a(a) takes the form 

*«) = (29) 

with 

a = A + 2 D q (30) 

For the present model Ao = 0. In fig. (4) we plot a(a) for different maps and compare with 
the Gaussian form given in eq.(29). For a constant h, non-Gaussian character increases 

6 



with increasing flatness of the map. But, as in the case of diffusion coefficient, increasing 
h has more influence in producing non-Gaussian character of the PDF. 

To conclude, analysis of the PDF with fluctuation spectrum brings out that intermit- 
tency and non-Gaussian character of the PDF increases with increasing peak height and 
flatness of the map, height exercising more effect. This is important when one selects 
models for describing experiments relating to diffusion. For example, systems exhibiting 
chaotic motion similar to Brownian motion should have a Gaussian distribution. Maps 
with linear segments having constant slope and minimum peak height are useful in cases 
like this. Maps with greater height with peaks becoming more flat will be best suited 
in describing diffusion systems showing intermittency in time. With flatness becoming 
maximum, diffusion coefficient behaves like 4r, h being the peak height. Corresponding 
limiting forms for characteristic function and fluctuation spectrum are also obtained. The 
limiting form of fluctuation spectrum is quite different from the Gaussian form following 
from central limit theorem. 

We are attempting a generalization of maps with fractional heights given in ref.(7) 
along similar lines. This work will be reported elsewhere. 

We are grateful to University Grants Commission, India, for providing financial assis- 
tance through DSA and COSIST schemes. 
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FIGURE CAPTIONS 



Fig. 1. Generalised Piecewise Linear (GPL) map with a variable shape with h — 1. 
On both axes, units are arbitrary. 

Fig. 2. Reduced map of GPL map in fig.(l). On both axes, units are arbitrary. 

Fig. 3. Variation of Diffusion coefficient D with r for h = 2. On both axes, units 
are arbitrary. 

Fig. 4. Fluctuation Spectrum a(a) versus a for different cases. Solid lines represent 
actual a (a) while dotted lines give corresponding Gaussian forms, (a) h — 1, m = 
mi = 7 (b) h — 1, m = 100, mi = 4.1237 (c) h = 2, m = mi = m 2 = 11. 
On both axes, units are arbitrary. 
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